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I .  Introduction 


This  report  describes  position  estimation  procedures  that 

t 

are  based  on  models  which  relate  positional  uncertainty  to 
various  measurement  and  estimation  errors.  The  procedures  were 
developed  to  be  used  in  analyses  that  relate  the  effect  of  posi¬ 
tional  uncertainty  on  tactical  performance  through  such  factors 
as  weapon  accuracy  as  well  as  to  be  used  operationally. 

In  the  models,  positions  are  on  a  plane  surface  (flat  earth 
model) .  Because  of  this  condition,  the  models  are  not  intended 
for  use  in  situations  where  the  earth's  figure  is  significant. 

In  addition,  positional  errors  are  determined  by  independent 
normally  distributed  random  variables  with  known  means,  variances 
and  covariances.  The  support  for  this  condition,  other  than  its 
mathematical  convenience,  is  that  it  has  been  used  by  others, 
for  example,  see  Reference  1.  To  use  the  models,  one  is  required 
to  specify  the  means,  variances  and  covariances  of  the  error 
random  variables.  For  these  models,  this  can  be  done  by  speci¬ 
fying  the  systematic  errors  (biases)  and  the  error  magnitudes 
(standard  deviations) . 

The  procedure  that  is  described  in  Section  II  relates  position 
estimates  that  are  based  on  bearings  on  or  from  stations  to 
station  position  uncertainty.  One  application  for  the  model  is 
the  analysis  of  the  effect  of  sonobuoy  position  uncertainty  on 
position  estimates  determined  with  passive  directional  sonobuoys . 

The  procedure  that  is  described  in  Section  III  relates  posi¬ 
tion  estimates  that  are  based  on  lines  of  position  to  line  of 


position  uncertainty.  In  addition  to  estimating  the  uncertainty 
in  position  estimates  based  on  celestial  observations,  the  pro¬ 
cedure  could  be  used  to  determine  error  ellipses  for  LORAN  fixes 
if  the  standard  deviation  values  required  by  the  procedure  could 
be  obtained. 

The  procedure  that  is  described  in  Section  IV  can  be  used 
to  combine  position  estimates  from  various  sources.  The  proce¬ 
dure  which  is  based  on  conditions  that  should  not  be  too  restric¬ 
tive  in  most  cases  provides  both  a  composite  position  estimate 
and  error  ellipse. 


II.  Station  Position  Uncertainty  and  Position  Estimates 

The  procedure  that  is  described  in  this  section  relates 
position  estimates  that  are  based  on  bearings  on  or  from  stations 
to  station  position  uncertainty.  The  procedure  is  based  on  a 
model  that  is  an  extension  of  one  that  is  described  on  Reference 
2.  The  model  is  defined  as  follows:  Each  station  position  error 
is  determined  by  an  independent  bivariate  normal  distribution 
with  a  zero  mean  vector  and  a  known  covariance  matrix.  Observed 
bearing  lines  on  or  from  a  station  are  parallel  to  true  bearing 
lines.  The  distance  of  each  observed  bearing  line  from  its  true 
bearing  line  is  determined  by  an  independent  normal  distribution 
with  a  zero  mean  and  a  standard  deviation  o. 

Because  a  station's  position  error  is  determined  by  a  bivar¬ 
iate  normal  distribution,  the  perpendicular  distance  s  between  a 
line  at  the  assumed  location  of  an  observed  bearing  line  and  the 
observed  bearing  line  is  determined  by  a  normal  distribution. 

The  relation  between  the  bivariate  normal  distribution  that 
describes  the  station's  position  uncertainty  is  indicated  in 
Figure  1.  In  the  figure,  the  positive  y-axis  direction  is  north, 
the  positive  x-axis  direction  is  east,  and  the  origin  of  the 
coordinate  system  is  at  the  assumed  station  position.  The  x'v'- 
coordinate  system  is  oriented  so  that  the  positive  x'axis  is 
coincident  with  the  major  axes  of  the  elliptical  contours  of 
the  bivariate  normal  distribution  that  determines  the  station 
position  error  and  so  that  the  bearing  6  of  the  positive  x'-axis 


Figure  1.  The  geometry  associated  with  the  determination  of  os, 
the  standard  deviation  of  the  normal  distribution 
that  determines  the  distance  between  the  assumed 
location  of  an  observed  bearing  line  and  the  observed 
bearing  line.  The  assumed  (mean)  position  of  the 
station  is  at  the  origin.  The  ellipse  represents  a 
contour  on  the  probability  density  surface  of  the 
bivariate  normal  distribution  that  determines  the 
station* s  position. 


satisfies  the  condition:  0°  <  5  <  180°.  The  angle  0  is  an 

observed  bearing.  The  x"y "-coordinate  system  is  oriented  so 
that  the  positive  y"-axis  is  in  the  direction  of  the  observed 
bearing  line.  As  a  consequence  of  these  relationships, 

os2  =  ax,2  sin2  ( 6-  0)  +  ay(2  cos2(6-e) 

2  2  . 
where  ax,  and  a^,  are  the  elements  of  the  station  position 

error  covariance  matrix  relative  to  the  x ' y ' -coordinate  system. 

A  procedure  for  determining  position  estimates  that  are 
based  on  bearings  on  or  from  stations  at  known  positions  is 
described  in  Appendix  1  of  Reference  3.  The  procedure  is  based 
on  a  model  that  is  equivalent  to  one  that  is  described  in  Ref¬ 
erence  2.  The  procedure  that  is  described  in  this  section  is 
based  on  a  model  that  is  an  extension  of  it. 

The  procedure  in  Appendix  1  of  Reference  3  is  based  on  a 
model  in  which  a  station's  bearing  error  is  determined  by  a 
normal  distribution  with  zero  mean  (bias)  and  standard  devia¬ 
tion  e.  The  bearing  error  is  related  to  the  distance  on  a  cir¬ 
cular  arc  between  a  station's  true  bearing  line  and  observed 
bearing  line.  The  arc  is  on  the  circle  with  its  center  at  the 
station  that  passes  through  an  initial  estimate  of  an  object's 
position.  This  distance  is  determined  by  a  normal  distribution 
with  mean  zero  and  standard  deviation  a  =  re  where  r  is  the 
range  of  the  initial  estimate  from  the  station  and  the  standard 
deviation  e  is  measured  in  radians.  In  the  model,  arc  distance 


is  approximated  using  a  first  order  approximation  which  in 
effect  replaces  the  circle  with  its  tangent  line  at  the  initial 
estimate.  As  a  consequence,  the  distance  u  on  the  tangent  line 
between  the  observed  bearing  line  and  the  true  bearing  line  is 
determined  by  a  normal  random  variable  with  mean  zero  and  stan¬ 
dard  deviation  a.  The  distance  u  can  be  expressed  in  terms  of 
w,  the  distance  on  the  tangent  line  between  the  observed  bearing 
line  and  the  initial  estimate  that  is  also  determined  by  a  normal 
random  variable  with  standard  deviation  a,  and  v,  the  distance 
on  the  tangent  line  between  the  true  bearing  line  and  the  initial 
estimate.  And,  as  shown  in  Appendix  1  of  Reference  3,  this 
distance  can  be  expressed  in  terms  of  the  unknown  coordinates 
of  the  object's  position. 

The  effect  of  station  position  uncertainty  is  accounted  for 
by  the  distance  s  between  the  observed  bearing  line  and  the 
assumed  observed  bearing  line.  In  the  model  s  is  determined  by 
a  normal  distribution  with  mean  zero  and  standard  deviation  a 

s 

as  given  above.  This  approximation  is  consistent  with  the  first 
order  approximation  of  arc  distance.  As  a  consequence  of  these 
two  approximations,  all  of  the  bearing  lines  are  replaced  by 
lines  parallel  to  the  line  joining  the  initial  estimate's 
position  and  the  station's  assumed  position,  both  of  which  are 
known  positions.  The  geometry  involved  in  shown  in  Fiqure  2. 

The  modified  relationships  resulting  from  the  introduction  of 
station  position  uncertainty  are  shown  in  Figure  3.  The 


assumed  location 
of  the  observed  bearing  line 


Figure  3.  The  quantities  u,  v  and  w  correspond  to  the  quantities 
u,  v  and  w  in  Appendix  1  of  Reference  3.  The  auxiliary 
quantities  u'  and  w'  are  defined  in  the  figure.  The 
replacement  of  by  +  og2  in  the  procedure  in 
Appendix  1  of  Reference  3  is  justified  by  noting  that 
with  this  substitution  u'  and  w'  are  equivalent  to  u 
and  w. 


3 


important  thing  to  note  is  that  u'  =  u  -  s  is  determined  by  a 

normal  distribution  with  mean  zero  and  standard  deviation 
2  2  k 

[a  +  o_  j  and  that  otherwise  U'  is  equivalent  to  U  with 

S 

respect  to  the  procedure  in  Appendix  1  of  Reference  3.  As  a 

con sequence  of  this,  the  procedure  can  be  extended  to  include 

2  2  h 

station  uncertainty  by  replacing  o  by  [o  +  o  ]  1  where  ever 
it  is  used.  In  this  case,  o  =  re  where  r  is  the  range  of  an 
initial  estimate  of  an  object's  position  from  the  assumed  station 
position  and  e  is  the  bearing  error  (standard  deviation)  in 
radians  of  the  bearings  associated  with  the  station. 


X  v.v.\ 


III.  Line  of  Position  Uncertainty  and  Position  Estimates 

The  procedure  that  is  described  in  this  section  relates 
position  estimates  that  are  based  on  lines  of  position  to  line 
of  position  uncertainty.  The  procedure  is  based  on  a  model  that 
is  defined  as  follows:  Lines  of  position  are  straight  lines. 
Observed  lines  of  position  are  parallel  to  true  lines  of  posi¬ 
tion.  The  distance  of  an  observed  line  of  position  from  a  true 
line  of  position  is  determined  by  an  independent  normally  dis¬ 
tributed  random  variable  with  known  mean  and  standard  deviation. 
Lines  of  position  are  specified  in  terms  of  a  rectangular  coor¬ 
dinate  system  with  the  origin  at  a  reference  point  as  shown  in 
Figure  4.  For  celestial  navigation,  an  appropriate  choice  for 
the  reference  point  would  be  the  assumed  position.  Since  bear¬ 
ing  lines  are  lines  of  position,  this  model  differs  from  the 
model  that  is  described  in  Section  II  only  in  terminology. 
However,  operationally  the  use  of  the  model  that  is  described 
in  this  section  differs  in  the  way  the  standard  deviation  of  the 
distance  of  the  line  from  a  true  line  is  determined.  The  stan¬ 
dard  deviation  associated  with  each  line  of  position  must  be 
specified.  If  this  can  be  done,  the  procedure  can  be  used. 

As  an  example,  suppose  that  the  values  are  o-^  for  the  first 
line  of  position  and  02  for  the  second  where  >  a 2  and 
that  the  lines  of  position  intersect  at  a  90°  angle.  In  this 
example,  the  minimum  area  confidence  (probability)  region  is 
an  ellipse  that  is  centered  on  the  estimated  position.  And  the 


ft. 
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Figure  4.  The  geometry  when  only  two  lines  of  position  are  used 
with  the  procedure  in  Appendix  1  of  Reference  3.  The  lines  of 
position  correspond  to  bearing  lines  for  observed  bearinqs  o,  and 
0£  with  respect  to  the  procedure.  The  procedure  determines  the 
size  and  the  orientation  of  elliptical  confidence  regions  of  con¬ 
fidence  p,  the  location  of  the  initial  estimate  and  the  location 
of  the  final  estimate.  The  locations  of  the  estimates  correspond 
when  only  two  lines  are  used.  If  the  lines  were  obtained  from 
sextant  observations,  the  assumed  position  should  be  the  reference 
position  in  which  case  the  lines  of  position  would  be  determined 
by  the  azimuth  anges  and  a ^  and  the  distances  pj  and  p2* 


estimated  position  is  at  the  intersection  of  the  bearing  lines. 

.  .  2 

For  a  confidence  (probability)  of  containment  of  1  -  exp(-k  /2) , 

the  major  axis  of  the  ellipse  is  coincident  with  the  first  bear¬ 
ing  line  and  it  is  of  length  2ko^,  the  minor  axis  of  the  ellipse 

is  coincident  with  the  second  bearing  line  and  it  is  of  length 

2 

2k0j  and  the  area  of  the  ellipse  is  irk  * 


IV.  A  Composite  Position  Estimate 


The  procedure  that  is  described  in  this  section  is  for 
combining  position  estimates  for  an  object  from  independent 
sources.  It  is  based  on  the  following  model:  The  rectangular 
coordinates  of  each  position  estimate  are  determined  by  an  inde¬ 
pendent  bivariate  normal  distribution  whose  covariance  matrix 
is  known  but  whose  mean  vector  is  not  known.  The  components  of 
the  mean  vector  for  each  of  the  distributions  are  the  unknown 
coordinates  x  and  y  of  the  object.  This  model  implies  that  the 
natural  logarithm  of  the  likelihood  function  for  a  set  of  n 
estimates  can  be  expressed  as  follows: 

n  ,  -1 

log  L  =  K  -  1/2  Z  (x.  -  x)  E.  (x.  -  x) 

i-1  “  1  “1 

where  K  is  a  constant,  x^  is  an  estimate  vector  with  components  x^ 

and  y^,  x  is  the  common  mean  vector  with  components  x  and  y, 

the  unknown  coordinates  of  the  object,  and  E^  is  the  covariance 

2  2 

matrix  with  elements  and  o .  The  maximum  likeli¬ 

hood  estimates  x  and  y  of  the  unknown  coordinates  x  and  y  are 
the  solutions  of  the  two  simultaneous  linear  equations  deter¬ 
mined  by 


3 (log  L) 
- 


y=y 


0  and 


3 (log  L) 

3y  x=x 

y=y 


0 


■A 
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The  equations  can  be  written  as 


A  x  +  B  y  =  D 
B  x  +  C  y  =  E 


and  their  solution  as 


x  = 


CD  -  BE 


AC  -  B 


-  AE  -  BD 

y  =  - 7 


AC  -  B 


where  A  =  Ea^,  B  =  Eb^,  C  =  Ec^,  D  =  E  (a^jL  +  b^y^) , 


E  -  Ubixl  +  c^l,  ax  =  ci?  /di;  b±  =  -  o±i~/di,  ct  =  o-  /di 


2  2  2 

d.  ■  o.-  -  o.0C,  and  all  of  the  sums  are  for  i  from  1  to  n 

x  lx  ly  xxy 


Since  they  are  linear  combinations  of  the  estimates  x^  and 

y^,  the  estimates  x  and  y  are  determined  by  a  bivariate  normal 

distribution.  Consequently,  all  that  is  required  to  determine 

this  distribution  is  its  mean  vector  with  components  y-  and  y~ 

2  2 

and  its  covariance  matrix  with  elements  o~  ,  o~  and  a  ~  ~ .  The 

x  y  xy 

mean  vector  is  determined  by 


y-  -  E{(CD  -  BE) /(AC  -  B2)}  =  x 

and 

y-  =E{(AE  -  BD)  /  (AC  -  B2)}  =  y 


And,  the  covariance  matrix  is  determined  by 


0-2  =  E{ (CD  -  BE)  -  E (CD  -  BE)}2/ (AC  -  B2)2 

=  (C2 (F+I+2H)  -  2CB(H+K+M)  +  B2 (G+J+2K)  }/ (AC  -  B2)2 


=  E{ (AE  -  BD)  -  E(AE  -  BD) }  /(AC  -  B  ) 


■  (A2(G+J+2N)  -  2AB (H+K+M)  +  B2 (F+I+2L) }/(AC  -  B2)2 


=  E ( [ (CD-BE) -E (CD-BE) ] * [ (AE-BD) -E (AE-BD) ] }/ (AC  -  B  ) 
xy 

*  { (AC+B2) (H+K+M) -CB(F+I+2L) -BA (G+J+2N) }/(AC  -  B2)2 

where  F  -  Ea/o^2  ,  G  =  Eb/o^2,  H  -  Ea^o^2  ,  I  =  Eb/o^2 

J  =  Ec/ck-2,  K  =•  Eb^o-2,  L  ■=  Ea^o.-.  M  -  E  (a^th/)  o^- 

and  N  =  Ib.c.o.~~  where  all  of  the  sums  are  for  i  =  1  to  n. 
xxx  xy 

By  using  arguments  given  in  Appendix  1  of  Reference  3,  one 

can  show  that  the  axes  of  the  elliptical  confidence  regions 

associated  with  x  and  y  are  coincident  with  an  x'y '-coordinate 

system  where  the  transformation  from  the  xy-coordinate  system 

to  this  system  is  the  coordinate  axes  rotation  through  the 

2  2 

angle  y  defined  by  tan  2y  =  -  o*  ) .  For  a  confidence 

xy  y  x 

p,  from  Reference  3,  the  minimum  area  confidence  region  is  an 

2 

ellipse  with  semi-axes  k  and  k  o*  •  and  area  ttJi 

x  y  #  x  y 

where  k  =  [-2  ln(l-p)]**  ,  , 

2  2  2  2  2 

o- ,  =0-  cos  y  -  2o~~  cos  y  sin  y  +  o~  sin  y, 

x  x  xy  y 

2  2  2  2  2 

o~,  =  o~  sin  y  +  20^  cos  y  sin  y  +  a~  cos  y 

y'  x  xy  y 

and  the  center  of  the  ellipse  is  at  the  point  (x,y).  In  this 


coordinate  system,  o~,*,  =  0 


The  above  equations  can  be  used  to  specify  a  composite  posi 

tion  estimate  in  terms  of  the  location,  orientation  and  size  of 

an  elliptical  confidence  region  which  is  generally  the  form  in 

which  position  estimates  of  the  kind  that  are  being  considered 

2  2 

here  are  specified.  But,  since  values  of  ,  0,7  and  ar.r. 

x  y  xy 

are  required  for  each  of  the  n  estimates  that  is  being  combined 
a  way  is  needed  for  determining  these  values  given  the  orienta¬ 
tion  and  size  of  an  elliptical  confidence  region.  A  procedure 
to  do  this  when  the  orientation  is  given  in  terms  of  the  direc¬ 
tion  6  of  the  major  axis  and  the  size  is  given  in  terms  of  the 
lengths  SMJ  and  SMI  of  the  semi-major  and  semi-minor  axes  and 
the  confidence  p  is  described  next. 

By  using  an  xy-coordinate  system  in  which  the  positive 
y-axis  direction  is  north  and  the  positive  x-axis  direction  is 
east  and  with  the  convention  0°  <  6  <  180°,  the  dependence  of 
the  value  of  the  rotation  angle  y  and  of  the  order  relation 

between  o~,  and  o~,  on  the  value  of  the  major  axis  direction  6 
x  y 

is  indicated  by  the  following  table: 

0°  <  6  <  45°:  y  -  6  and  o~,  >  a*, 

“  y  * 

45°  <  6  <  135°:  y  =  6  -  90°  and  o~,  >  o~, 

“  x  y 

135°  <  6  <  180°:  y  ■  6  -  180°  and  a-,  >  o~, 

With  an  order  relation  and  a  value  for  p,  values  for  o*,  and 

op,  can  be  determined  with  values  for  SMJ  and  SMI.  With  values 

2  2 

for  0£,  ,  o'),  and  y,  values  for  0*  ,  o^  and  can  be  found 

x  y  x  y  xy 

from  the  following  equations: 


°xy  =  ^ay '  “  °x*  *  sin  Y  cos  Y 

which  can  be  obtained  by  inverting  the  equations  above  for  , 

,  2 
and  o- , 

As  an  example,  suppose  the  data  in  the  following  table 
represent  three  independent  position  estimates: 


A 

X 

A 

y 

6 

SMJ 

SMI 

k 

1st 

-3.7 

18.1 

59° 

36 

20 

2 

2nd 

11.8 

• 

CD 

105° 

37 

11 

2 

3rd 

0 

0 

146° 

45 

23 

2 

Here, 

distances 

are  in 

nautical 

miles  and  p  = 

.86  in  each 

case 

For  this  example  with  values  in  square  nautical  miles: 


X> 

K) 

II 

264. 

.58, 

°iy2  = 

159. 

.42 

and 

o  = 

lxy 

98.89 

(1 

CM 

<>i 

CM 

321. 

.35, 

°2y2  = 

51. 

.15 

and 

°2xy  = 

-78 

2 

A  ss 

3x 

249. 

K) 

O 

2 

o,A  = 
3y 

389. 

.30 

and 

°3xy  = 

-173.4 

These  values  give  the  following  composite  estimate: 
x  =  -2.46  nautical  miles  and  y  =  12.26  nautical  miles. 


For  this  case,  a~  =  16.92  nautical  miles,  o~  =  6.21  nautical 

x  y 

miles  and  =  -55.1  square  nautical  miles.  And,  for  k  =  2: 
xy 

SMJ  =  33.85  nautical  miles,  SMI  =  12.42  nautical  miles  and 
6  =  102°.  The  composite  confidence  region  and  its  three  com¬ 
ponent  confidence  regions  are  shown  in  Figure  5. 

As  a  second  example,  suppose  each  position  estimate  is  deter 
mined  by  a  circular  normal  distribution.  Then 

am1,  o./v)  =  0  for  i  =  1  to  n.  In  this  case,  the  composite 
ixy 

estimate  is: 

x  =  (^i/oi)/(El/o.)  ,  y  =  (Eyi/oi)/(Zl/oi)  , 

a~2  »  n/(El/a.)2  ,  a-2  =  n/ (El/a.) 2  and  a**  =  0. 
x  i  y  i  xy 

In  this  example,  since  x  and  y  are  determined  by  a  circular 
normal  distribution,  the  minimum  area  confidence  regions  are 
circles  and  orientation  is  not  an  issue. 


North 


East 

x 


Figure  5.  The  ellipses  represent  the  confidence  ellipses  of 
the  first  example.  The  position  estimates  are  at  the  center 
of  the  ellipses.  The  numbers  indicate  the  order  of  the  esti¬ 
mate  in  the  table  on  Page  17. 
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